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ABSTRACT 

Using the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, 
we have obtained high-resolution ultraviolet observations of GD 362 and PG 1225- 
079, two helium-dominated, externally-polluted white dwarfs. We determined or 
placed useful upper limits on the abundances of two key volatile elements, car- 
bon and sulfur, in both stars; we also constrained the zinc abundance in PG 
1225-079. In combination with previous optical data, we find strong evidence 
that each of these two white dwarfs has accreted a parent body that has evolved 
beyond primitive nebular condensation. The planetesimal accreted onto GD 362 
had a bulk composition roughly similar to that of a mesosiderite meteorite based 
on a reduced chi-squared comparison with solar system objects; however, addi- 
tional material is required to fully reproduce the observed mid-infrared spectrum 
for GD 362. No single meteorite can reproduce the unique abundance pattern 
observed in PG 1225-079; the best fit model requires a blend of ureilite and 
mesosiderite material. From a compiled sample of 9 well-studied polluted white 
dwarfs, we find evidence for both primitive planetesimals, which are a direct 
product from nebular condensation, as well as beyond-primitive planetesimals, 
whose final compositions were mainly determined by post-nebular processing. 

Subject headings: planetary systems stars: abundances white dwarfs 

1. INTRODUCTION 

Planetesimals are building blocks of planets and their formation is a key step towards 
planet formation. How do planetesimals form? What determines their bulk composition? 
To answer these questions, we start by examining our own solar system. 
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The overall configuration of the solar system is that volatile-depleted, dry rocky objects 
are ubiquitous relatively close to the Sun while volatile-rich, icy objects are found beyond 
the snow line. This correlation between the volatile fraction and heliocentric distance can 
be explained by primitive nebular condensation: refractory elements condensed closer to 
the Sun while volatile elements can only be incorporated into the planetesimals where the 
temperature is low enough. Many solar system objects have experienced some additional 
processing that changed their initial compositions. For example, it has been argued that a 
collision between a large asteroid and proto- Mercury s tripped off most of Mercury's silicate 
mantle, leaving it ~70% iron by mass (IBenz et al. I ll988l). Also, the "late veneer" has delivered 
a large amount of water and volatiles onto Earth (jChybal Il990l ). Post-nebular processing, 
such as collisions, melting and differentiation, is important in redistributing the elements 
among solar system objects. 

Currently, the best way to measure the elemental compositions of planetesimals in the 
sol ar system is from meteor ites, which are fragments from collisions among asteroids. Follow- 
ing |O^NeilL&.£almg (120081 ). we classify all meteorites into two categories in this paper, (i) 
"Chondritic" is used to refer to chondrites, which are a direct product of nebular processing. 
Objects in this category are described as "primitive" planetesimals. (ii) "Non-chondritic" 
objects consist of achondrites, stoney-iron meteorites and iron meteorites. Examples of their 
parent bodies include the Moon, Mars or asteroids that have experienced various amounts 
of post-nebular processing. Planetesimals in this category are considered to be "beyond- 
primitive" . 

What about planetesimal formation in extrasolar planetary systems? High-resolution, 
high-sensitivity spectroscopic observations of externally-polluted white dwarfs ar e a powerfu l 
tool for determining the bulk elemental compositions of extrasolar planetesimals ( lJurall2013l ). 
Calculations show that minor planets can survive the red giant stage of a star and p ersist 
into the white dwarf phase with most of their internal water and volatiles intact (jjura 
20081 : 1 Jura fc Xull2010l ). Orbital perturbations from one or multiple planets can cause these 
planetesimals to stray into the tidal radius o f the white dwarf and get tidally disrupted 
( iDebes fc Sigurdsson 



2002; Bonsor et al. 2011: Debes et al. 



dust disk that enaits naostly in the infrar ed (jJura 



2003 



2012[). sometimes producing a 



Kilic et al.ll2006l : Ivon Hippel et al. 



20071 : iFarihi et al.ll2009l : IXu &: Jural |2 0121 ). Eventually, all this planetary debris is accreted 



onto the central white dwarf and pollutes its otherwise pure hydrogen or helium atmosphere. 

The first com prehensive abundance rn easurement of an externally-polluted white dwarf 
was performed by IZuckerman et al.l ( 20071 ). who identified 15 elements heavier than helium 
in the atmosph ere of CD 362 , including Mg, Si and Fe, which are often called the "com- 
mon elements" ( lLarimerlll988l ). The disrupted object had a minimum mass ~10^^ g, which 
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is comparable to that of a massive solar system asteroid. Three years later, the abun- 
dances of eight heavy elements were determined in the atmosphere of GD 40, including all 
the major rock-forming elements - O, Mg, Si and Fe ( iKlein et al.ll2010l ). Now there are 
more high- r esolut io n optical sp e ctrosc o pic studies of exte r nally- p olluted white dwarfs 



many more mgn- r esoiut io n optical sp e ctrosc o pic studies oi exte r naiiy- p oiiuted wnite dwaris 
[e.g.. iKlein et all fcoilh: iMelis et all fcoilhJzuckerman et all fcoilh : IParihi et all fcoilj ): 



Dufour et all (120121 ): IVennes et all feoiol . 



20111 )]. 



However, optical spectroscopy of externally-polluted white dwarfs typically does not 
enable sensitive detection of highly-volatile elements, such as carbon, nitrogen and sulfur, 
which are key to understanding the thermal history of the system. Ultraviolet spectroscopy 
is complimentary to optical observations in determination of volatile abundances. 

To-date, there are four white dwarfs with both published high-resolution optical and 
ultraviolet mea^wementfl we are beginning to accumulate an atlas of the compositions 
of extrasolar planetesimals. To zeroth order, we find that they are strikingly similar to 
meteorites in the solar system: (i) O, Mg, Si and Fe are always dominant and their sum 
is more than 85% of the accreted mass; (2) volatile elements, especially C, are typically 
depleted by more than a factor of 10 compared to solar abundance^. 

In this paper, we report ultraviolet spectroscopic observations of GD 362 and PG 1225- 
079 with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope (HST), 
complimentary to previous optical studies from the Keck High Resolution Echelle Spectrom- 
eter (HIRES) dZuckerman et ahlboO?! : iKlein et al.lboilh . PG 122 5-079 has been o bserved 
with the low- resolution International Ultraviolet Explorer (lUE) (jWolff et al.l |2002[ ) : there 
is no previous ultraviolet spectroscopy for GD 362. The rest of the paper is organized as 
follows. Data reduction is summarized in section 2 and atmospheric abundance determi- 
nations are reported in section 3. In section 4, we used a reduced chi-squared analysis to 
look for solar system analogs to the accreted parent bodies. The formation mechanisms 
of extrasolar planetesimals are assessed in section 5 and conclusions are given in section 6. 
In Appendix A, we report the Herschel Photodetecting Array Camera and Spectrometer 
(PACS) observation of GD 362. In Appendix B, we extend the reduced chi-squared analysis 
to two additional externally-polluted helium white dwarfs with both high-resolution optical 
and ultraviolet observations. 



m 



^The four wh i te dwarfs are: GD 61 (IDesharnais et al.l booa IParihi et all boill): GD 40. G241-6 



Klein et al.ll20ld '2011: 'Zuckerman et al."20ld: ljura et al.l l2012h and WD 1929+012 ( Vennes et aUboiol 



2011i : .Melis et al...2011 : Gansickc et al..,20121 



^Very recently, iKoester et al.l (j2012l) reported several white dwarfs with solar carbon-to-silicon ratio. 
However, the source of this pollution is unclear and more analysis is forthcoming. 
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OBSERVATIONS AND DATA REDUCTION 



GD 362 and PG 1225-079 were observed during HST/ COS Cycle 18 under program 
12290. These two white dwarfs are too cool to be observed ef f ective ly w ith the G130M 



gratin g centering around 1300 A, as was employed by lJura et al.l ( |2012| ) and lCansicke et al. 



( 120121 ) for other hotter white dwarfs. Instead, the G185M grating was used with a central 
wavelength of 1921 A and wavelength coverage of 1800 - 1840 A, 1903 - 1940 A and 2008 - 
2044 A. The spectral resolution was ~18,000. Total exposure times were 7411 and 1805 sec 
for GD 362 and PG 1225-079, respectively. 

The raw data were processed using the standard pipeline CALCOS 2.13.6. The fluxes 



at 2030 A are 2.9 x 10"^^ erg s"^ cm'^ A"^ and 1.5 x 10"^^ erg s"^ cm'^ A"^ for GD 362 
and PG 1225-079, respectively, in approximate agreement with broadband NUV fluxes from 
the GALEX satellite. The signal-to-noise ratio (SNR) in the original un-smoothed spectrum 
was 6 for PG 1225-079 and 4 for GD 362. 



Following previous data reduction procedures ( iKlein et al.ll2010l l2011t I Jura et al.ll2012l ). 
for PG 1225-079, equivalent widths (EWs) of each spectral line were measured in the un- 
smoothed spectra by fitting a Voigt profile with three different nearby continuum intervals 
in IRAF. The EW uncertainty is calculated by adding the standard deviation of the three 
EWs and the average uncertainty from the profile fitting in quadrature. The EW upper limit 
is obtained by artificially inserting a spectral line with different abundance into the model 
and comparing with the data. We adopt a different method to measure the EW for C I 
1930.9 A in GD 362, as described in section 3.1. The measured values are listed in Tables 
[Hand [2] for GD 362 and PG 1225-079, respectively. The average Doppler shift relative to 
the Sun for PG 1225-079 is 42 ± 13 km s~^, in essen tial agreement with the value 49 ± 3 
km s~^ derived from optical studies (IKlein et al.l 1201 ll ). The large velocity dispersion in the 
ultraviolet is due to the low SNR of the spectrum and the ~ 15 km s~^ uncertainty of COS 
(COS Instrument Handbook). For GD 362, we marginally detected C I 1930.9 A and it has 
a Doppler shift of 48 km s~^, in agreement with 49.3 ± 1.0 km s~^ from the optical study 
( IZuckerman et al.l 120071 ). 



3. ATMOSPHERIC ABUNDANCE DETERMINATIONS 



Because we are most interested in the abundance of an element relative to other heavy 
elements and t hese ratios are no t strongly dependent upon the stellar temperature and 
surface gravity ( iKlein et al.l 1201 ll ) , we only adopt one set o f stel lar parameters as listed in 
Table [3] and compute the model spectra following iKoesterl ( 120101 ). Atomic data are mostly 
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Table 1: Measured Equivalent Widths and Abundance Determinations for GD 362 



Ion 


A 

(A) 


(eV) 


EW 
(mA) 


log n(Z)/n(He) 


C 1 


1930.905 


1.26 


560 tlTs " 


-6.70 ± 0.30 


S 1 
S I 

S 


1807.311 
1820.341 




0.049 


<900 
<710 


< -6.70 

< -6.40 

< -6.70 



" This is measured from the model spectra, as described in section 3.1. 



taken from the Vienna Atomic Line Database (iKupka et al.lll999l ). The computed model 
atmosphere spectra were convolved with the COS NUV line spread functioiif|. The abundance 
of each element was derived by comparing the EW of each spectral line with the value 
derived from the model atmosphere, as shown in Figures [TH5] and Tables [1] and O The final 
abundances, combining ultraviolet with optical observations, are given in Tables H] and [5] 
for GD 362 and PG 1225-079, respectively. Our results mostly agree with previous reports 
but have a higher accuracy. For PG 1225-079, we newly derive the abundances of carbon 
and silicon and have tentative detections of sulfur and zinc. The magnesium abundance is 
updated while the iron abundance agrees with previous optical results. Because the data 
are noisier for GD 362, we are only able to crudely constrain the abundance of carbon and 
sulfur. 



3.1. Carbon 



There is only one useful carbon line in the observed wavelength interval, C I 1930.9 
A, as shown in Figures [1] and [2l Because it arises from an excited level, it cannot be 
contaminated by interstellar absorption. However, this line can be blended with Mn II 
1931.4 A. For tunately, accurate Mn a bundances have b een determined for both stars from 
optical data ( IZuckerman et al.l 120071 : iKlein et al.l 1201 ll ) and the predicted EW for Mn II 
1931.4 A is less than 50 mA in the model spectrum. Considering the measured EW of this 
feature is more than 500 mA for both stars (see Tables [1] and [2]) , we conclude that the line 



http://www.stsci.edu/hst/cos/perforniaiice/spectral_resolution/nuv_modeLlsf 
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Table 2: Measured Equivalent Widths and Abundance Determinations for PG 1225-079 



Ion 


A 




EW 


log n(Z)/n(He) 




(A) 


(eV) 


(mA) 




C I 


1930.905 


1.26 


1600 ± 200 


-7.80 ± 0.10 


S I 


1807.311 





< 170 


< -9.50 


S I 


1820.341 


0.049 


< 150 


< -9.30 


S 








< -9.50 


Mg I 


2026.477'^ 





288 ± 100^ 


< -7.60 


O * TT 

bi 11 


1808.013 





936 ± 109 


''7 A A 1 n "1 M 

-7.44 ± 0.10 


Si II 


1816.928 


0.04 


1232 ± 145 


-7.46 ± 0.10 


Si 








-7.45 ± 0.10 


Fe II 


1925.987 


2.52 


192 ± 72 


-7.62 ± 0.28 


Fe II 


2011.347 


2.58 


309 ± 97 


-7.35 ± 0.24 


Fe II 


2019.429 


1.96 


211 ± 67 


-7.56 ± 0.24 


Fe II 


2021.402 


1.67 


181 ± 68 


-7.71 ± 0.27 


Fe II 


2033.061 


2.03 


311 ± 67 


-7.24 ± 0.17 


Fe II 


2041.345 


1.964 


215 ± 50 


-7.24 ± 0.18 


Fe 








-7.45 ± 0.23 


Zn II 


2026.136 





288 ± 47^ 


< -11.30 



" The atomic parameters for this line are taken from iKelleher fc Podobedoval (120081 ). 

^ Mg I 2026.5 A and Zn II 2026.1 A are blended and the reported EW is for the entire 

feature. 
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Table 3: Adopted Stellar Properties 



star 


M, 


T 


log g 


D 


log M,,,/M/ 


Ref 




(Me) 


(K) 


(cm^ s^"*^) 


(pc) 






GD 362 


0.72 


10,540 


8.24 


51 


-6.71 


(1) (2) 


PG 1225-079 


0.58 


10,800 


8.00 


26 


-5.02 


(3) (4) 



" Newly-derived mass of the conve ctive zone (see section 4). 

References , f 11 iKilic et al.l f|2008h : (2) IZuckerman et all fl2007l ): (3) iKlein et al.l J201lh : (4) 



Farihi et all fl2005[ ). 



is dominated by C I 1930.9 A. For PG 1225-079, our derived carbon ab undancjjl [C]/[He 



-7.80 ± 0.10 agrees with the lUE upper limit of -7.5 flWolff et al.l 120021 ). 



For GD 362, the largest uncertainty is from the low SNR of the data; the measured 
continuum flux is (3.1 ± 1.0) x 10~^^ erg s~^ cm~^ A~^. It is hard to measure the EW of C 
I 1930.9 A directly from the noisy data. Instead, we computed model spectra with different 
carbon abundance to match the observed spectrum. In Figure [11 we present three best-fit 
models with [C]/[He] = - 6.4, [C]/[He] = -6.7, [C]/[He] = -7.0 and a continuum flux at 4.1 
X 10"-^^ erg s~^ cm"^ A"\ 3.1 x 10"-^^ erg s"-^ cm"^ A~^, 2.1 x 10~^^ erg s"^ cm"^ A~^, 
respectively. The final abundance is [C]/[IIe] = -6.7 ± 0.3 and the EW reported in Tabled] 
is measured from the model spectra. 



3.2. Sulfur 



There are two useful sulfur lines, S 1 1807.3 A and S 1 1820.3 A. However, at best, we have 
only a tentative detection of sulfur in each star. S I 1807.3 A, the stronger line, is adjacent to 
Si II 1808.0 A. Fortunately, for GD 362, the silicon abundance is determined from previous 
optical data (IZuckerman et al.l 120071 ): for PG 1225-079, other ultraviolet lines can be used 
to derive the silicon abundance (see section 3.4). The data and model atmosphere spectra 
for GD 362 and PG 1225-079 are presented in Figures [3] and HI respectively. Considering the 
apparent match between the model and data for both S I lines, tentative sulfur abundances 
of -6.7 for GD 362 and -9.5 for PG 1225-079 can be assigned. Conservatively, these results 



■^Here, log n(X)/n(Y) is abbreviated as [X]/[Y]. 
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Table 4: Atmospheric Abundances for GD 362 



rj 


log n(/j/n(Hej 




b 

et 


M[LiY 






(10^ yr) 


(g s ) 


TT 

H 


-1.14 ± U.IU 










-0.7U ± U.3U 


2 


1 


2.5 X 10 


1\,T 


^ AAA 

< -4.14 


2 


2 


< 9.0 X 10^ 


(J 


< -5.14 


2 


2 


< 1.1 X 10 


AT 

IMa 


-7.79 ± 0.20 


2 


2 


3.7 X 10 


Mg 


-5.98 ± 0.25 


2 


2 


2.5 X 10 


A 1 
Al 


-D.40 ± 0.20 


1 


6 


1.5 X 10 


bi 


r G A 1 n on 

-5.84 ± O.oO 


1 


2 


7.2 X 10 


b 


< -6.70" 





79 


< 1.7 X 10 


Ca 


-6.24 ± 0.10 





99 


5.1 X 10 


Sc 


-10.19 ± 0.30 





93 


6.8 X 10^ 


Ti 


-7.95 ± 0.10 





94 


1.2 X 10^ 


V 


-8.74 ± 0.30 





95 


2.1 X 10^ 


Cr 


-7.41 ± 0.10 


1 





4.3 X 10^ 


Mn 


-7.47 ± 0.10 


1 





4.0 X 10^ 


Fe 


-5.65 ± 0.10 


1 


1 


2.5 X 10*^ 


Co 


-8.50 ± 0.40 





99 


4.1 X 10^ 


Ni 


-7.07 ± 0.15 


1 





1.1 X 10^ 


Cu 


-9.20 ± 0.40 





83 


1.1 X 10^ 


Sr 


-10.42 ± 0.30 





56 


1.3 X 10^ 


Total 








4.4 X 10^ 



* New measurements from this paper. The rest are from IZuckerman et al.l (120071 ) but we 
reference abundances r elative to He, the domin ant element in GD 362's atmosphere, rather 
than H, as presented in lZuckerman et al.l (120071 ). Consequently, there is a possible systematic 
offset up to 0.1 dex in all entries derived from that paper. 

The final abundance of an element combining optical and ultraviolet data. 
^ Newly-derived settling times in the convective zon e (see section 4) ; they are typically a 
factor of 2-3 longer than previously- derived values in iKoesterl (120091 ). 

Accretion rates calculated from Equation (1). 

The equality sign corresponds to the red model fit shown in figures. 
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Table 5: Atmospheric Abundances for PG 1225-079 



z 


log n(Z)/n(He) 


tset 








(10^ yr) 


(g s-i) 


H 


-4.05 ± 0.10 






C* 


-7.80 ± 0.10 


5.5 


3.1 X 10^ 





< -5.54 


4.5 


< 9.1x 10^ 


Na 


< -8.26 


4.4 


< 2.6 X 10^ 


Mg* 


-7.50 ± 0.20 


4.8 


1.4 X 10^ 


Al 


< -7.84 


3.6 


< 9.5 X 10^ 


Si* 


-7.45 ± 0.10 


3.0 


3.0 X 10^ 


S* 


< -9.50 


1.7 


< 5.2 X 10^ 


Ca 


-8.06 ± 0.03 


1.9 


1.6 X 10^ 


Sc 


-11.29 ± 0.07 


1.8 


1.1 X 10^ 


Ti 


-9.45 ± 0.02 


1.8 


8.3 X 10^ 


V 


-10.41 ± 0.10 


1.8 


9.6 X 10^ 


Cr 


-9.27 ± 0.06 


1.9 


1.3 X 10^ 


Mn 


-9.79 ± 0.14 


2.0 


4.0 X 10^ 


Fe 


-7.42 ± 0.07 


2.1 


9.0 X 10^ 


Ni 


-8.76 ± 0.14 


2.3 


4.0 X 10*^ 


Zn* 


< -11.30 


2.2 


< 1.3 X 10^ 


Sr 


< -11.65 


1.2 


< 1.4 X 10^ 


Total 






1.6 X 10* 



* New results from this paper. The rest are from lKlein et al, 
Notes. The columns are defined the same as Table HI 



(1201 If ). 
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Fig. 1. — HST/COS spectrum of GD 362. The black line is the data smoothed with a 3 pixel 
boxcar. The green, red, blue line represents the computed model spectrum with [C]/[He] = 
-6.4, -6.7, -7.0, respectively, placed at a different continuum level; the abundances of other 
elements are from Table |H The adopted carbon abundance is -6.7 ± 0.3. The red labels 
represent lines that are used for abundance determinations. Wavelength is presented in the 
star's reference frame in vacuum. 




1922 1924 1926 1928 1930 1932 1934 1936 1938 1940 



A(A) 



Fig. 2. — Similar to Figure [T] except for PG 1225-079 with abundances from Table El The 
data were smoothed with a 5 pixel boxcar. 



are upper limits. 



3.3. Magnesium and Zinc 



In PG 1225-079, Mg I 2026.4 A and Zn II 2026.1 A are heavily blended. As shown 
in Figure our best fit model which matches the measured EW of the absorption feature 
requires [Mg]/[He] = -7.6 and [Zn]/[IIe] = -11.3. These values are individually taken as 
upper limits due to the ble nding. How'ever, the reported magnesium abundance is -7.27 ± 
0.06 from the optical data ( iKlein et al.ll201l|). which i s large ly based on three Mg lines but 
the detections for two lines are only 2a. IWolff et al.l ( 120021 ) reported [Mg]/[He] to be -7.6 
=b 0.6 from the lUE data. Averaging these measurements, our final magnesium abundance 
is -7.50 ± 0.20. Because of the blending, the zinc abundance is only an upper limit. This 
provides the first stringent constraint on zinc in an extrasolar planetesimal. 



3.4. Silicon 



In PG 1225-079, we measured two silicon lines. Si II 1808.0 A and Si II 1816.9 A, as 
shown in Figure |H Si II 1808.8 A arises from the ground state and the photospheric line 
can be distorted by interstellar absorption. However, its measured EW is only 87 ± 11 mA 
i n C Qph; st ar at a distance of 112 pc with a large amount of foreground interstellar gas 
( iMortonI 119751 ). Considering PG 1225-079 is only 26 pc away, it has much less interstellar 
absorption. The measured EW is 936 ± 109 mA and we conclude that Si II 1808.0 A is 
largely photospheric and essentially free from interstellar absorption. The shape of Si II 
1808.0 A in the model does not quite fit the data; but the measured EW of the data, which 
is key in the abundance determination, has a good agreement with that in the model. Using 
these two Si II lines, we derive a final silicon abundance of -7.45 ± 0.10, in agreern ent with, 
but much better than the reported l UE abundance of - 7.5 ± 0.5 ( IWolff et al.l 120021 ) and the 
previous optical upper limit of -7.27 (IKlein et al.l 1201 if ). 



3.5. Iron 



In the COS data for PG 1225-079, there are six Fe II lines with EWs larger than 100 
mA. Four of them are shown in Figures [2] and [5l We derived an iron abundance of -7.45 
± 0.23, in good agreement of the optical value of -7.42 ± 0.07, which is based on 28 high- 
SNR iron lines (IKlein et al.l 1201 ll ). Because the ultraviolet data are noisier, we adopt the 
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1825 



Fig. 3. — HST/COS spectrum of GD 362. All notations are the same as Figure [T] and the 
data are smoothed by a 5 pixel boxcar. S I 1807.3 A and 1820.3 A are used for constraining 
the sulfur abundance. 




1815 

A(A) 



1825 



Fig. 4.— Similar to Figure [3] except for PG 1225-079. Si II 1808.0 A and 1816.9 A lines are 
used for determining the silicon abundance. 
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2025 

A(A) 



2035 



Fig. 5. — HST/COS spectrum of PG 1225-079 and the data are smoothed by a 5 pixel 
boxcar. All notations are the same as Figure [H Zn II 2026.1 A and Mg I 2026.5 A are used 
for determining the Zn and Mg abundances, respectively. The green line presents the best 
fit model that matches the overall EW of the absorption feature at 2026 A ([Mg]/[He] = 
-7.6, [Zn]/[He] = -11.3); the red line is the adopted model combining with the optical results 
([Mg]/[He] = -7.5, [Zn]/[He] = -11.3). Fe II 2019.4 A, 2021.4 A, 2033.0 A and 1926.0 A in 
Figure [2] are used for determining the iron abundance. 



optically-derived iron abundance. 



COMPARISON WITH SOLAR SYSTEM OBJECTS 



Combined with previous data, we now have determined the abundances of 16 elements 
heavier than helium in the atmosphere of GD 362 and 11 heavy elements in PG 1225-079. 
However, the measured composition need not be identical to the composition of the accreted 
planetesimal because different elements gravitationally settle at different rates in a white 
dwarf atmosphere. Thr ee major phases are proposed for a single accretion event: build-up, 
steady-state and decay (IDupuis et al.lll993l : lKoesterll2009l ). 



Because an infrared excess is found for GD 362 and PG 1225-079 ( Becklin et al.l 12005 



Kilic et al.ll2005t iFarihi et al.ll2010l ). the accretion should be either in the build-up o r steady 
state phase. The timescale for build-up stage is comparable to the settling times ( iKoester 
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20091); it is ~ 10^ yr, for GD 362 and PG 1225-079 (see Tables 1 and [5]) . The rest of the 
disk-host stage should all be under the steady-state approximation. The dust disk lifetime 
has been under intensive studies for a few year s but t he values are still very uncertain 
including 10^ yr jP arihi et al.ll2009H R afikovll2011b[) , 10^ y r jRafikovll2011al : lGirven et al.ll2012 



Farihi et al. 20121 ) and up to 10^ yr ( Barber et al. 2012). The true disk lif etime might have 



a rang e but it is likely to be longer than the settling times. Furthermore, IZuckerman et al. 



(I2OIOI ) suggested that steady-state approximation is the dominant situation for white dwarf 
accretion event based on a study of helium dominated stars; the settling times are only 
0.1% of their cooling times but 30% of them show atmospheric pollution. GD 362 and PG 
1225-079 are more likely to be under the steady-state approximation and that is the main 
focus of this paper. 

In the steady-state model, the observed concentration of an element is dependent on 
the time it takes to sink out of the convective envelope. To derive the theoretical settling 
times and obtain an improved understanding of the uncertainties, we formulated several 
numerical experiments with the code for the envelope structure and corrected two errors 
found in our previous calculations of diffusion timescales . In the course of ch anging the 
equations descri bing element di f fusion from the version in iPaquette et al.l (119861 ) (Equation 
4) to the one in iPelletier et al.l ( 119861 ) (Equation 5), which is more accurate in the case of 
electron degeneracy, one of us (D.K.) discovered an error in the former paper. A factor of 
pi/3 jg jnissing in the second alternative of Equation 21, which we had not noticed before. 
A rederivation of all our equations uncovered another error in our implementation of the 
contribution of thermal diffusion. These errors have only a very small effect in stars with 
relatively shallow convection zones, like the hydrogen-dominated white dwarfs. However, for 
helium-dominated white dwarfs with T < 15,000 K and a deep convection zone, the diffusion 
timescales can be slower by factors 2-3 relat ive to ou r earlfe. calculationfl The accretion 
rate M{Zi) of an element Z is calculated as (IKoesteii 120091 ) 



M{Zi 



:ii 



where Mc^^ is the mass of the convective envelope. X(Zj) is the mass fraction of the element 
Zj relative to the dominant element in the atmosphere, either hydrogen or helium; tsei(Zj) is 
the settling time. A longer settling time corresponds to a lower diffusion flux. Fortunately, 
the relative timescales for different elements, which are important for the determination of 
the abundances in the accreted matter, change much less. 



^Updated diffusion timescales can be obtained at 



http://www.astrophysik.uni-kiel.de/~koester/astrophysics/ 
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For GD 362 and PG 1225-079, compared to previously published values, the settling 
times listed in Tables H] and |5] typically increase by factors of 2-3 while the mass of the 
convective zone is 0.13 dex smaller for GD 362 and 0.05 dex larger for PG 1225-079 (Table 
in]). These corrections lead to smaller total accretion rates by a factor of 3 for both stars. 

The next step is to compare the composition of the accreted parent body with those of 
solar system objects. We choose the summed mass of all the major elements as the normal- 
ization factor so that the analysis is independent of the chemical property and abundance 
uncertainty of each individual element. However, one complication is that no oxygen lines 
are detected in either GD 362 or PG 1225-079 due to their low photospheric temperatures 
relative to other helium-dominated white dwarfs; only upper limits were obtained for this 
major element. Therefore, our approach is to compare the mass fraction of an element rela- 
tive to the summed mass of the common elements Mg, Si and Fe. For solar sy stem objects. 



we inc lude 80 representative and well-analyzed meteorite samp les mostly fro r n iNittler et al. 
(2004). We also include the bul k compositio n of Earth from AUegre et al. ( 200ll ) and an 



updated carbon abundance from iMartyl (120121 ). For our purpose. Earth appears to be chon- 
dritic and its bulk composition approaches CV chondrites even though Earth has experienced 
some post-nebular processing, such as differentiation and collisions. 



4.1. GD 362: Accretion from a Mesosiderite Analog? 

In Figure [6l we compare the abundances of all 18 elements, including upper limits, of the 
accreted material in GD 362 with CI chondrites, which are the most primitive material in the 
solar system. The composition of CI chondrites is almost identical to the solar photosphere, 
with the exception of depletion of volatile elements C, N as well as H and noble gases. The 
parent body accr e ted o nto GD 362 looks nothing like a CI chondrite, as first pointed out in 
Zuckerman et al.l (120071 ). For the volatile elements, the mass fraction of C and S are depleted 
by at least a factor of 7 and 3, respectively, relative to CI chondrites; refractory elements, 
such as V, Ca, Ti and Al, are all enhanced. 

Though oxygen is n ot detected i n GD 362, its stringent upper limit can still provide 
useful insights. Following iKlein et al.l (120101 ). we can calculate the required number of oxygen 



atoms to form oxides 'Zip(^z)Og(z) as 



(2) 



Hydrogen is excluded here because GD 362 has an enormous amount and it might not be 
associated with the parent body or bodies currently in its atmosphere (see Appendix A). 
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Fig. 6. — Mass fractions of heavy elements in GD 362 with respect to the summed mass of 
sihcon, iron and magnesium from Table HI The abundances are normalized to those of CI 
chondrites. The elements are ordered by decreasing volatility. The filled stars represent the 
steady-state approximation while the open stars are in the build-up stage, la error bars 
and arrows for upper limits are plotted only for the steady-state approximation for clarity. 
We also plot the avera ge values for five me sosiderites (Emergy, Barea, Patwar, Dyarrl Island 
and ALH 77219) from lNittler et al.l ( 120041 ): la deviations are shown as blue error bars. The 
error bars are smaller than the symbol when they are not visible in the figure. There are no 
reported whole rock fractions of Cu, V, Sr and Sc in mesosiderites. 



Under the steady-state approximation, [0]/[He] = -5.07 is required to form MgO, AI2O3, 
Si02 and CaO; this value is comparable to the observed oxygen upper limit of -5.14. However, 
Fe is the most abundant heavy element in the atmosphere of GD 362 and there is insufficient 
oxygen to tie it up in either FeO or Fe203. Thus, most, if not all the iron in the parent body 
i s in metallic form, which is very different from CI chondrites where most iron is in oxides 
dNittler et al]l2004l ). 



O'Neill fc Palmd (120081) suggested that [Mn]/[Na] can be used as an indicator of post- 



nebular processing. For example, [Mn]/[Na] is -0.79 for all chondrites as well as the solar 
photosphere while non-chondritic objects have a much higher value. Interestingly, [Mn]/[Na] 
i s 0.65 ± 0.22 for GD 362, which is larger than -0.01 for Mars and 0.32 for the Moon 
(jO'Neill fc Palmd I2OO8I ). This suggests that the planetesimal accreted ont o GD 362 is likely 



to be non-chondritic and have experienced some post-nebular processing. IZuckerman et al. 



(I2OO7I ) compared the [Na]/[Ca] ratio in GD 362 with solar system objects and reached a 
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similar conclusion; the accreted planetesimal was non-chonridtic. The only other polluted 
white d warf with both Mn and Na detections is WD J0738+1835 wherein [Mn]/[Na]= -0.54 



± 0.19 ( iDufour et al.ll2012[ ): this agrees with the chondritic value within the uncertainties. 



To find the best solar system analog to the parent body accreted onto GD 362, we 
calculated a reduced chi-squared value for each object in our sample {Xred)y defined as: 

where N is the total number of elements considered in the analysis. Mu,(^(Zj) and Mmtri'Z'i) 
represent the mass fraction of an element Zj relative to the summed mass of Mg, Si and Fe in 
the extrasolar planetesimal and solar system objects, respectively. a^dCZii) is the propagated 
uncertainty in mass fraction. 

For GD 362, we calculated Xred ^ heavy elements, C, Na, Mg, Al, Si, Ca, Ti, Cr, 
Mn, Fe and Ni, which have detections both in GD 362 and the meteorite sampl^. The 
results are shown in Figure [7] for both steady-state and build-up approximations. There 
is no qualitative difference between these two models and mesosiderites provide the best 
fit considering all 11 elements. In particular, the mesosiderite ALH 77219 can match the 
overall abundance pattern to 95% confidence level. As shown in Figure [6l the abundance of 
individual elements agrees within 2a between mesosiderites and the planetesimal accreted 
onto GD 362. 

Mesosiderites are a rare type of stoney-iron meteorite with equal amounts of silicates 
and metallic iron and nickel. One mystery about mesosiderties is that the Si-rich crust and 
Fe, Ni-rich core materials are abundant but the olivine Mg-rich mantle seems to be missing. 
One model for the formation of mesosiderites is that a 200-400 km diameter asteroid with 
a molten core was nearly catastroph ically disrupted by a 50-150 km diameter projectile at 



4.42-4.52 Gyr ago (IScott et al.ll200ll ). The collision mixed the target's molten core with its 
crustal material but excluded the large and hot mantle fragments. The planetesimal accreted 
onto GD 362 may have been formed in a similar way. 

While mesosiderites may be a prototype for the accreted planetesimal onto G D 362, there 



are th ree major hurdles for this hypothesis to overcome. First, in the model of IScott et al. 



(l200ll ). only half of the original mass of a 200-400 km diameter asteroid was maintained 



after the collision and the final product only contains about 10% mesosiderite-like material 



^For a couple of meteorites with no reported carbon abundance, we compute the Xred ^'^^ ^^^^ other 10 
elements. 
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Fig. 7. — This figure shows x^^^ values defined in Equation (3), which compares the abun- 
dances of C, Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe and Ni in the accreted material in GD 
362 with 80 meteorites and bulk Earth. The x-axis denotes Xred calculated for the compo- 
sition under the steady-state approximation while the y-axis is for the build-up phase. The 
dashed lines represent 95% confidence level. There are a couple of meteorites with a large 
Xred and they are not shown in the current scale. There is no qualitative difference between 
the steady-state and build-up models. Mesosiderites, particuarly ALH 77219, provide the 
best fit to the abundance pattern observed in GD 362. Notes: "Chondritic" materials in- 
clude: (i) carbonaceous chondrites: CI, CK, CM, CO, CR and CV; (ii) ordinary chondrites: 
H, L and LL; (iii) enstatite chondrites: EH and EL; (iv) R chondrite. "Non-chondritic" 
materials consist of achondrites and stoney-iron meteorites. No iron meteorite is included 
in our analysis because their bulk composition is dominated by Fe and Ni, with very few 
trace elements reported. Achondrites include: (i) primitive achondrites: acapulcoites, lo- 
dranites, winonaites and ureilites; (ii) asteroidal meteorites: angrites, aubrites, brachinites 
and howardite-eucrite-diogenite; (iii) Martian meteorites: shergottie, Nakhlites and Chas- 
signites; (iv) Lunar meteorites. For stoney-iron r neteorites, we includ e (i) mesosiderites; (ii) 
pallasites. Most of the meteorite data are from iNittler et al.l ( 120041 ) and the compositions 
for s ome Martian nieteor ites are from iMcSweeru (119851 ). The bulk composition of E arth is 
from lAUegre et al.l (120011 ) and the carbon abundance is from a mor e recent study of iMartv 
( 120121 ) , which is a factor of 3 lower than the lower limit reported in lAUegre et al.l ( 12001! ) . 
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by mass. This is equivalent to a 75-150 km diameter object. Mesosiderites that fall on Earth 
are only small fragments and the 180 kg NWA 2924 is among the largest (Meteorite Bulletin 
Databas^). However, the parent body accreted onto GD 362 has a minimum mass of 2.7 
X 10^^ g, ~260 km in diameter for an assumed density of 3 g cm~'^. It is unclear whether 
the same kind of collision can produce a mesosiderite parent body this big. Second, the 
mass fraction of hydrogen in mesosiderites is less than 0.2%; it cannot explain how there is 
5 X 10^^ g hydrogen in the atmosphere of GD 362. Possibly, hydrogen was accreted during 
earlier events and it has been atop the atmosphere ever since (see Appendix A for more 
discussion). Third, GD 362 is currently accreting from its circumstellar disk and the disk 
material should also resemble the composition of mesosiderites. However, the shape of the 
mid-infrared spectrum for mesosiderite, whi ch is dominated by a sharp peak at 9.13 fim and 
several other bands at 10.6 /xm and 11.3 /im (IMorlok et al.ll2012l). cannot fully account for the 
broad 10 fim silicate emission feature observed for GD 362 (I Jura et al.ll2007l ). This does not 
completely exclude the mesosiderite hypothesis but emission from some additional material 
is required to fully reproduce the observed infrared spectrum for GD 362. Mesosiderites 
are a good candidate for the parent body accreted onto GD 362 but there are remaining 
unresolved issues. 



4.2. PG 1225-079: Accretion from a Planetesimal with No Single Solar 

System Analog 

In Figure [8], we show a comparison of the mass fractions of 16 elements, including upper 
limits between PG 1225-079 and CI chondrites. Though the carbon abundance is approaching 
the chondritic value, the accreted planetesimal differs a lot from CI chondrites; the mass 
fraction of S is depleted by at least a factor of 40 while Zn is depleted by at least a factor of 
8. In contrast, refractories, such as V, Ca, Ti and Sc are all enhanced. The overall pattern of 
relatively high carbon abundance and enhanced mass fractions of refractory elements does 
not follow a single condensation sequence and post-nebular processing is required. 

As shown in Figure [H PG 1225-079 has a [C]/[S] value that is no smaller than the solar 
ratio, which is very different from other polluted white dwarfs and meteorites. Carbon and 
sulfur are among the most volatile elements that we can measure and their 50% condensation 



temperatures are 40 K and 655 K, respectively (jLoddersll2003l ). Most of the meteorites as 
well as polluted white dwarfs have a [C]/[S] ratio lower than the solar value, which can 
be explained by condensation at a temperature between 40 and 665 K though this is not 



^http://www.lpi.usra.edu/nieteor/ 
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(b) 

Fig. 8. — (a) Similar to Figure |6] except for PG 1225-079 with the abundances in Tabled 
Ureihtes can match the high carbon low sulfur patt ern in PG 1225-079 but fail with the 
other elements. The abundances for ureilites are from lWarren et al.l (120061 ). (b) The best-fit 
model for the composition in PG 1225-069 under the steady-state approximation is 30% by 
mass ureilite North Haig and 70% mesosiderite Dyarrl Island. There are no reported bulk 
compositions of Zn, V or Sc for the North Haig and Dyarrl Island meteorites. 
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necessarily true for all of them. The only solar system analog to PG 1225-079 with similar 
high carbon, low sulfur pattern is ureilites, a type of primitive achondrites. Ureilites are the 
second largest achondrite group and it is suggested that its high carbon abundance is derived 
from a carbon-r ich parent body, but the exact formation mechanism is not well understood 
( lGoodriclj|l992l ). However, as can be seen in Figure |8]^a), ureilites fail to match the overall 
composition of the parent body accreted onto PG 1225-079. 

We performed a xled analysis between solar system objects and the accretedplanetesimal 
in PG 1225-079, comparing 9 elements, C, Mg, Si, Ca, Ti, Cr, Mn, Fe and NO- The result 
is shown in Figure [101 There is no single solar system object that can match all nine 
elements; the closest is carbonaceous chondrite. Regardless, as shown in Figure El the 
accreted abundance in PG 1225-079 is not at all identical to CI chondrites. 



The infrared excess around PG 1225-079 corresponds to ~500 K dust ( iFarihi et al. 



20101); so far, only two white dwarfs are known to have s uch cool dust. The other 28 known 



disk- host stars all have ~1000 K dust (jXu fc Jural 120121 ). One hypothesis is that the inner 



Farihi et al. 


2010; 


Jura 


2008) 



of two planetesimals, rather than one single parent body. After testing different combinations 
of the 80 meteorites in our database, the best fit model to the steady-state approximation 
consists of 30% ureilite North Haig and 70% mesosiderite Dyarrl Island by mass. This blend 
is also marked in Figure [TUl Detailed abundance comparison is shown in Figure IHI^b); the 
abundances of S, Mn and Ca do not agree as well as the other elements but are all within 
2a. A possible scenario is that one extrasolar ureilite (mesosiderite) analog first got tidally 
disrupted and more recently, another mesosiderite (ureilite) analog impacted the disk and 
was blended with the previous material. 



5. ASSESSING THE FORMATION MECHANISMS OF EXTRASOLAR 

PLANETESIMALS 

Having established that the parent bodies accreted onto CD 362 and PG 1225-079 
are beyond primitive, we now extend our analysis to other extrasolar planetesimals. We 
are most interested in understanding the formation mechanisms of extra solar planetesimals 



and whether these are dominated by nebular or post-nebular processing. lJura fc Xul (120131 ) 



suggested collisional rearrangement is important in determining the final composition of 



^Similar to the case of GD 362, for the meteorites with no reported carbon abundance, we only calculated 
xled the other 8 elements. 
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Fig. 9. — Mass fraction of S and C over the sum of Fe, Mg and Si for solar system objects as 
well as polluted white dwarfs with positive detection of carbon or sulfur or both. Due to the 
presence of an infrared excess except for G241-6 (see discussion in Appendix B), all white 
dwarfs are plotted under the steady-state approximation . la uncertainties are plotted. The 
red dashed line denotes constant solar [C]/[S]. Most meteorites have a lower [C]/[S] than 
the solar value. However, [C]/[S] in PG 1225-079 is no smaller than solar and the closest 
solar system analog is ureilites. Refere nces: WD 1 9 29-F0 12, PG 0843+517, WD 1226+110: 



Gansicke et all (j2012|); GD 40. G2 41-6: Jjura et all (120121 ): GD 362 and PG 1225-079: this 



paper; Solar abundance: iLodderd (120031 ): for solar system objects, the references are listed 
in the caption of Figure [71 
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Fig. 10. — Similar to Figure [7] except for PG 1225-079 comparing 9 elements, C, Mg, Si, 
Ca, Ti, Cr, Mn, Fe and Ni. We see that no solar system object comes close to the overall 
abundance pattern in PG 1225-079. The black triangle represents 30% by mass ureilite and 
70% mesosiderite as shown in Figure [H](b), which is the best match for the composition of 
PG 1225-079 in the steady-state approximation. 



extrasolar planetesimals based on the scatter in [Mg]/[Ca] ratios in 60 externally-polluted 
white dwarfs. Here, we compile a sample of well-studied externally-polluted white dwarfs 
with abundance determinations of at least 9 elements. There are 9 stars in total, as listed 
in Table inland now we assess the formation mechanism for individual objects. 



GD 40: As discussed in lJura et al.l ( 120121 ) and Appendix B, the overall abundance pat- 
tern in GD 40 matches with carbonaceous chondrites and bulk Earth. Nebular condensation 
is sufficient to explain its observed composition. 



WD J0738+1835: iDufour et al.l ( 120121 ) found that there is a correlation between the 
abundance of an element and its condensation temperature: refractory elements are de- 
pleted while volatile elements are enhanced compared to bulk Earth. This indicates that 
the accreted planetesimal might be formed in a low temperature environment under nebular 
condensation. 



Table 6: Summary of Planetesimal Formation Mechanisms in 9 Well-Studied White Dwarfs 



star 


Dom. 


Dust 


Volatile 


Intermediate 


Refractory 


Process 


Ref 


GD 40 


He 


Y 


C,S:,0,Mn, P 


Cr,Si,Fe,Mg,Ni 


Ca,Ti,Al 


primitive 


1,2 


WD J0738+1835 


He 


Y 


0,Na,Mn 


Cr,Si,Fe,Mg,Co,Ni 


V,Ca,Ti,Al,Sc 


primitive 


3,4 


PG 0843+517 


H 


Y 


C,S,0,P 


Cr,Si,Fe,Mg,Ni 


Al 


b eyond-pr imit i ve ( ? ) 


5 


PG 1225-079 


He 


Y 


C,Mn 


Cr,Si,Fe,Mg,Ni 


V,Ca,Ti,Sc 


beyond-primitive 


6,7 


NLTT 43806 


H 


N 


Na 


Cr,Si,Fe,Mg,Ni 


Ca,Ti,Al 


beyond-primitive 


8 


GD 362 


He 


Y 


C,Na,Cu,Mn 


Cr,Si,Fe,Mg,Co,Ni 


V,Sr,Ca,Ti,Al,Sc 


beyond- primitive 


7,9 


WD 1929+012 


H 


Y 


C,S,0,Mn,P 


Cr,Si,Fe,Mg,Ni 


Ca,Al 


??? 


5,10,11,12 


G241-6 


He 


N 


S,0,Mn,P 


Cr,Si,Fe,Mg,Ni 


Ca,Ti 


primitive 


2,6,13 


HS 2253+8023 


He 


N 


0,Mn 


Cr,Si,Fe,Mg,Ni: 


Ca,Ti 


primitive 


6 



Note. This is a compiled sample of externally polluted white dwarfs with detections of at least 9 elements heavier 
than helium. Columns are defined as follows. "Dom" lists the dominant element in the atmosphere. "Dust" indicate s 
whether a star has an infrared excess ("Y") or not ("N"). Following the classification scheme in iLodderd (120031 ) . 
"Volatile" lists the detected volatil e elements, de fined as having a 50% condensation temperature lower than 1290 K 
in a solar-system composition gas (ILoddersI l2003l ): "Intermediate" lists the elements with a condensation temperature 
between 1290-1360 K -the same range as that of the common elements. Si, Fe and Mg; "Refractory" elements 
have a 50% condensation temperature higher than 1360 K. The elements are ordered with increasing condensation 
temperature. "Process" shows our proposed dominant mechanism that determines the final composition of the accreted 
extrasolar planetes imal (see se c tion 5 ). 
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PG 084-3+517: This star has the highest mass fraction of iron among all polluted white 
dwarfs. iGansicke et al.l (|2012[ ) found that all core elements, including Fe, Ni, S and Cr are 
enhanced relative to the values for bulk Earth while lithophile refractory Al is depleted. This 
star might be accreting from the core of a differentiated object. Nevertheless, considering 
the uncertainty for each element is at least 0.2 dex, the conclusion is still preliminary. 

PG 1225-079: As discussed in section 4.2, this star has a near chondritic carbon abun- 
dance but also enhanced mass fractions of refractory elements relative to CI chondrite; it 
cannot be formed solely under nebular processing. 

NLTT 43806: Com pared to chondrit i c valu es, the accreted planetesimal is depleted in 
Fe and enhanced in Al. IZuckerman et al.l ( 120111 ) found that the best fit model corresponds 



to "30% crust 70% upper mantle". With detections of 9 elements, evidence is strong that 
NLTT 43806 has accreted the outer layer of a differentiated parent body. 

GD 362: As discussed in section 4.1, mesosiderite is the best solar system analog to the 
accreted parent body and post-nebular processing is required. 



WD 1929+012: IGansicke et al.l (120121 ) showed that this star has a high iron content. 
However, the situation is perplexing in that different analyses yie l d diffe rent stellar parame- 
ters a nd atmospheric abundances. For e xample. both lMelis et al.l (1201 ll ) and IGansicke et al. 
J2012h derived that [Si]/[Fe] is -0.25 but IVennes et al.l J201oh found that [Si]/[Fe] is 0.19. No 
final conclusion can be drawn before resolving such discrepancies. 

G24I-6: This star is a near twin of GD 40 with a similar abundance pattern but without 
an infrared excess. One possible scenario is that G241-6 has accreted a planetesimal with 
a similar composition to GD 40 and now it is at the beginning of a decaying phase; all 
heavier elernents appear to be dep leted relative to GD 40 due to the ir short setthng times 
(IKlein et al.ll201ll : I Jura et al.ll2012l ). As discussed in I Jura et al.l (120121 ) and Appendix B, the 
overall abundances resemble those of chondrites and no post-nebular processing is required. 



HS 2253+8023: iKlein et al.l (I2OIII ) showed that the composition of its parent body 
agrees with bulk Earth, except for the enhanced calcium abundance. Nebular processing 
can produce the observed abundance pattern. 

As summarized in Table [6l at least 4 out of the 9 white dwarfs have accreted planetesi- 
mals that can be formed under nebular processing while post-nebular processing is required 
for another 3 of them. It should be noted that some objects that we identify as primitive 
might still have undergone some post-nebular processing. For example, GD 40 has accreted 
from a planetesimal that has a similar composition as bulk Earth, whose overall abundance 
pattern is chondritic. However, it is still possible that the parent body was differentiated; 
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when the entire object is accreted, the composition appears to be "chondritic" . We can only 
put an upper hmit on the number of objects formed under nebular condensation. 

From this sample of 9 stars, we see that post-nebular processing appears to play an 
important role in determining the final abundance of extrasolar planetesimals; beyond- 
primitive planetesimals might be as common as primitive planetesimals. In contrast, chon- 
drites comprise more than 90% of all meteorites found on Earth by number (Meteorite 
Bulletin Databas^. Possibly, extrasolar planetesimals around white dwarfs have violent 
evolutionary histories with more collisions. This difference is not surprising since dynamical 
rearrangement of planetary systems at white dwarfs is expected to increase the frequency of 
collisions and produce more beyond-primitive extrasolar planetesimals. 

So far, 19 elements heavier than helium, including C, S, O, Na, Cu, Mn, P, Cr, Si, Mg, 
Fe, Co, Ni, V, Sr, Ca, Ti, Al and Sc, have been detected in the atmospheres of polluted 
white dwarfs, as shown in Table |6l In terms of mass fracti on in the accreted p lanetesimal, 
the lowest limit is ~5 ppm, for Sc in WD J0738+1835 JPufour et all l2012h . Studying 
externally-polluted white dwarfs proves to be a very sensitive probe of the bulk compositions 
of extrasolar planetesimals. 



6. CONCLUSIONS 

We present HST/COS ultraviolet observations for GD 362 and PG 1225-079, two heavily 
polluted helium white dwarfs. In GD 362, the mass fractions of carbon and sulfur are depleted 
by at least a factor of 7 and 3 respectively, compared to CI chondrites. In PG 1225-079, a 
similar volatile depletion pattern is found: C by a factor of 2, S by at least a factor of 40 and 
Zn by at least a factor of 8. We provide good evidence for the presence of beyond-primitive 
extrasolar planetesimals: 

1. Mesosiderites provide a good match to the composition of the parent body accreted 
onto GD 362. However, there are several unresolved issues for this hypothesis, espe- 
cially the apparent difference between the mid-infrared spectrum of mesosiderites and 
the dust disk around GD 362. Additional material is required. 

2. No single meteorite can reproduce the abundance pattern in PG 1225-079. A blend of 
30% North Haig ureilite and 70% Dyarrl Island mesosiderite can provide a good fit to 
the overall composition. 



http://www.lpi.usra.edu/nieteor/ 
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3. Spectroscopic observations of externally-polluted white dwarfs enable sensitive mea- 
surement of the bulk compositions of extrasolar planetesimals, including 19 heavy 
elements down to a mass fraction of 5 ppm. Based on a sample of 9 well-studied white 
dwarfs, we find that post-nebular processing is as important as nebular condensation 
in determining the compositions of extrasolar planetesimals. 

Support for program 7^ 12290 was provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. This work also has been partly 
supported by NSF grants to UCLA to study polluted white dwarfs. 

APPENDIX 



A. The Herschel/PACS Observation of GD 362 



While hydrogen is detected in some helium-dominated white dwarfs (jVoss et al.l 120071 ). 
GD 362 has an anomalously large amount. The helium-to-hydrogen number ratio is 14 in 
its convecti ve zone, correspon ding to 5 x 10^^ g of hydrogen; this is lower than 7 x 10^^ g 
reported in lJura et al.l (120091 ) because the mass of the convective zone for GD 362 is 0.13 
dex lower in the updated calculation (Table EJ. The origin of the hydrogen is a mystery. 
Unlike heavy elements which have short settling times compared to the white dwarf cooling 
age, hydrogen never s i nks and can be ac cumulated over the entire cooling history of the star 
( iBergeron et al.ll201ll : lJura fc Xul 120121 ). If GD 362 has always been a helium-dominated 
white dwarf and all this hydrogen is from accretion of ti dally disrupted o biects, it can either 
be one Callisto-size object or ~100 Ceres-like asteroids ( jJura et al.ll2009l ). In the latter case, 
likely there would be many more asteroids orbiting the star and mutual collisions among 
them would generate a cloud of cold dust. 



We were awarded 1.1 hours of Herschel/PACS (jPoglitsch et al.ll2010l ) observation time 
to look for cold dust around GD 362. The "mini-scan map" mode was used to observe in 
"blue" (85-125 /im) and "red" (125-210 /im) bands simultaneously with a medium scan speed 
of 20" s^^ and a scan leg length of 4'. The scan map size is 345" x 374" and the repetition 
number is 25. Two different scan angles, 45 degrees and 135 degrees were used and the total 
integration time was 1200 sec. 

Data reduction was performed using HIPE (Herschel Interactive Processing Environ- 
ment) on a combined mosaic of level 2 products from pipeline SPG 7.1.0. The pixel scale is 
1" pixeP^ and 2" pixeP^ for the blue and red band, respectively. Correcting for its proper 



- 28 - 



motion, we expect GD 362 at a = 17:31:34.355, S = +37:05:18.331 on the date of the obser- 
vation. Because there is no detection, aperture photometry was performed at 25 locations 
within 5 pixels of the nominal position of GD 362. The aperture radius was 20" with a sky 
annulus between 61" and 70". The background intensity was estimated using the median 
sky estimation algorithm (Herschel Data Analysis Guidq^. Aperture correction factors are 
0.949 for blue and 0.897 for red (PACS Observer's Manua0). Based on the dispersion of 
the 25 measurements, 3a upper limits are 5.1 mJy for blue and 5.6 mJy for red. 

What does this imply about dust mass? G D 362 has shrunk in mass from 3 Mq on the 
main-sequence to its current mass of 0.72 M© (IKilic et al.l l2008l ). Consequently, asteroids 
initially at 3-5 AU are now orbiting at 13 -21 AU. Currently , GD 362 has a stellar temperature 
of 10,540 K and c ooling age ~ 0.9 Gyr (IFarihi et al.ll2009l ). Extrapolating from white dwarf 
cooling modeli^ ( Bergeron et al. 2011 ). for GD 362, its stellar temperature is lower than 



20,000 K for 90% of its cooling time. We approximate the stellar luminosity as a time- 
averaged luminosity of 0.01 L©. Poynting-Robertson drag was able to remove particles 
smaller than 20 /im at a distance of 15 AU for a grain density of 3 g cm~^. We therefore 
assume a dust particle radius of 20 /im in the putative asteroid belt orbiting GD 362. 

If the grains function as blackbodies with negligible albedo, then their temperature can 
be calculated as 



2D orb 



(Al) 



T*, R* are the stellar temperature and radius; Dorb is the orbital distance. The dust 
temperature is 14-11 K between 13-21 AU. 

The mass of the dust disk is 



XBAT) 



(A2) 



where D^, is the distance to GD 362, 51 pc (IKilic et al.ll2008l ) and x is the dust opacity. For a 
particle radius of 20 /im, x = 100 cm^ g^^ in the geometric optics limit. As shown in Figure 
[TTl the upper limit of dust mass is between 10^^ g and 10^^ g at 13-21 AU; this mass is at 



^'^http://herschel.esac.esa.int/hcss-doc-8.0/print/howtos/howtos.pdf 
^^http: / /herschel.esac.esa.int /Docs /PACS / pdf/pacs_oni.pdf 
^^http://www. astro. umontreal.ca/ bergeron/CoolingModels/ 
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least twice the hydrogen mass in the atmosphere of GD 362 and one or der of magnitude 
larger than the mass of solar system's asteroid belt (IKrasinsky et al.ll2002l ). The upper limit 
is not stringent enough to rule out the hypothesis that hydrogen in GD 362 is from accretion 
of multiple asteroids. So, the large hydrogen abundance in GD 362 remains an unsolved 
puzzle. 




15 20 
Dorb (AU) 



30 



Fig. 11. — Upper limit for the mass of cold dust around GD 362 derived from PACS blue 
and red data, as a function of orbital radius. Given the hydrogen mass of 5 x 10^^ g, the 
upper limit of dust mass 10^^-10^^ g at 13-21 AU is not stringent enough to rule out the 
accretion of multiple asteroids. 



B. Looking for Solar System Analogs to Extrasolar Planetesimals 



The Xred analysis has proven to be an effective way to look for solar system analogs 
to the compositions of extrasolar planetesimals. Two other helium-dominated white dwarfs 
have reported volatile and refractory abundances from high-resolution optical and ultraviolet 
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observations that are suitable for this kind of analysia^^l - GD 40 and G241-6. Updated 
settling times and accretion rates are listed in Table [7] while the mass of the convective zone 
stays the same. Since all the major elements are determined, we compare the mass fraction 
of an element relative to the sum of O, Mg, Si and Fe. 



Table 7: Updated Settling Times and Accretion Rates for GD 40 and G241-6 





i. a 


M(Z) 


GD40 


M(Z)g241-6 


z 


(10^ yr) 


(g s"^ 


) 


(g s-^) 


c 


1.1 


2.2 X 


10^ 


< 4.4 X 10^ 


N 


1.1 


< 2.6 


X 10^ 


< 2.1 X 10^ 





1.1 


4.5 X 


10^ 


4.3 X 10^ 


Mg 


1.2 


1.7 X 


108 


1.5 X 10^ 


Al 


1.2 


1.4 X 


10^ 


< 6.1 X 10^ 


Si 


1.0 


1.3 X 


10^ 


8.7 X 10^ 


P 


0.79 


1.1 X 


10^ 


4.7 X 10^ 


s 


0.64 


1.0 X 


10^: 


5.6 X 10^ 


CI 


0.51 


< 8.0 


X 10^ 


< 5.8 X 10^ 


Ca 


0.51 


1.3 X 


10^ 


5.1 X 10^ 


Ti 


0.49 


3.2 X 


10^ 


1.4 X 10^ 


Cr 


0.53 


6.4 X 


10^ 


4.5 X 10^ 


Mn 


0.53 


3.1 X 


10^ 


2.4 X 10^ 


Fe 


0.56 


4.4 X 


10« 


2.0 X 10* 


Ni 


0.61 


1.8 X 


10^ 


8.9 X 10^ 


Cu 


0.58 


< 1.8 


X 10^ 


< 1.8 X 10^ 


Ga 


0.50 


< 2.9 


X 10^ 


< 2.9 X 10^ 


Ge 


0.43 


< 1.4 


X 10^ 


< 1.4 X 10^ 


Total 




1.4 X 


10^ 


9.9 X 10* 



" This column is for both GD 40 and G241-6 because their atmospheric conditions are similar. 



The total accretion rate for GD 40 is a factor of 2 lower than the value derived in 



Klein et al.l (l2010l ). but the relative abundances change much less. The result of a xled 
analysis is presented in Figure [121 When including all 13 detected elements, both carbona- 
ceous chondrites and bulk Earth can match the composition to 95% confidence level for 



^^G D 61 also has high-resolution optical and ultraviolet observations (|Desharnais et al.ll2008t iFarihi et al. 
20111 ). However, with a total of 5 detected elements, it is hard to make a comparison using the xled analysis. 
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both steady-state and build-up approximations. The accreted planetesimal appears to be 
primitive a r id can be formed under nebular condensation, similar to what was concluded by 



Jura et al.l fl2012f ). 
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Fig. 12. — Similar to Figure [7] except for GD 40 comparing the mass fraction of 13 elements, 
including C, O, Mg, Al, Si, P, S, Ca, Ti, Cr, Mn, Fe and Ni, relative to the sum of Mg, Si, 
Fe and O. Both carbonaceous chondrites and bulk Earth are a good match to the parent 
body accreted onto GD 40. 

The newly-deriv ed total accretion rate for G241-6 is about a factor of 2 lower than 



previously reported (jZuckerman et al.ll2010l ). The non-detection of an infrared excess and 



(Xu & Jura 


2012; 


Klein et al. 


2011) 



for the xied analysis; in the latter case, we assume that accretion stopped 0.6 x 10 yr 
ago, approximately one settling time for Fe because its mass fraction is depleted by a factor 
of 2 relative to CI chon drites. The compos ition of the pa rent body is calculated following 
Zuckerman et al.l ( l201ll ) and Equation (5) in lKoesterl (l2009l ). A fuller exploration of diff erent 
time- varying models will be presented in the future in the spirit of lJura fc Xul (|2012| ). As 
shown in Figure [131 both carbonaceous chondrites and ordinary chondrites provide good 
matches to all 11 elements, including O, Mg, Si, P, S, Ca, Ti, Cr, Mn, Fe and Ni. However, the 
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carbon upper limit in G241-6, which is not included in th e x^^^ analysis, is at least one order 
of magnitude lower than most carbonaceous chondrites ( iJura et al.ll2012l ). Thus, ordinary 
chondrites are a more promising solar system analog to the parent body accreted onto G241-6 
and nebular condensation is sufficient to produce the observed abundance pattern. 



The xled analysis for GD 40 and G241-6 conffims the previous results (I Jura et al.ll2012l ): 



the accreted extrasolar planetesimals can be formed under nebular condensation and their 
compositions resemble primitive chondrites in the solar system. 
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Fig. 13. — Similar to Figure [T2] except for G241-6 comparing 11 elements - O, Mg, Si, P, S, 
Ca, Ti, Cr, Mn, Fe and Ni for steady-state versus decay phase when the accretion stopped 
0.6 X 10^ yr ago. Both carbonaceous chondrites and ordinary chondrites produce a good fit 
to the parent body accreted onto G241-6. Ordinary chondrites are a relatively better match 
because of the low carbon abundance, w hich is not consid ered in this xled pl*^^ because only 
an upper hmit was reported for G241-6 ( I Jura et al.ll2012[ ). 
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